Graphene on the carbon face of SiC: electronic structure modification by hydrogen 
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It has been shown that the first C layer on the SiC(0001)(2 x 2)c surface already exhibits graphene- 
like electronic structure, with linear tt bands near the Dirac point. Indeed, the (2 x 2)c reconstruc- 
tion, with a Si adatom and C restatom structure, efficiently passivates the SiC(OOOl) surface thanks 
to an adatom/restatom charge transfer mechanism. Here, we study the effects of interface mod- 
ifications on the graphene layer using density functional theory calculations. The modifications 
we consider are inspired from native interface defects observed by scanning tunneling microscopy. 
One H atom per 4x4 SiC cell (5x5 graphene cell) is introduced in order to saturate a restatom 
dangling bond and hinder the adatom/restatom charge transfer. As a consequence, the graphene 
layer is doped with electrons from the substrate and the interaction with the adatom states slightly 
increases. Native interface defects are therefore likely to play an important role in the doping mecha- 
nism on the C terminated SiC substrates. We also conclude that an efficient passivation of the C face 
of SiC by H requires a complete removal of the reconstruction. Otherwise, at variance with the Si 
terminated SiC substrates, the presence of H at the interface would increase the graphene/substrate 
interaction. 



PACS numbers: 81.05.ue, 31.15.A- 



I. INTRODUCTION 
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Graphene, a single layer of graphite, has two C atoms 
per unit cell, usually called A and B, that form a hon- 
eycomb lattice. Its low-energy electronic structure ex- 
hibits two TT bands, the it band is filled and the tt* band 
is empty. They touch at the Dirac point, at the cor- 
ner of the Brillouin zone. The corresponding energy is 
called the Dirac Energy {Ed). Importantly, the disper- 
sion is linear within ±0.5eV with respect to i?zj[H. The 
experimental discovery of graphene's unique electronic 
structure f2l, and high electronic mobility]!] announces 
a whole new era for physics as well as nanoelectronics[B, 
0] . Among the fabrication techniques like mechanical ex- 
foliation from a graphite crystal [Ti] or catalytic decom- 
position of hydrocarbons on transition metals graphi- 
tization of the polar faces of SiC has the advantage of 
providing large scale graphene sheets, directly on an in- 
sulating substrate^. However, since graphene is an ul- 
timately thin crystal, its environment has to be care- 
fully considered. Proving that point, the highest mobil- 
ity so far has been achieved for suspended graphene 
Hence, for supported graphene, it is of particular impor- 
tance to investigate the interface structure and to study 
how the substrate impacts the electronic properties of the 
graphene layer. 

For graphene grown on the Si terminated substrates 
(SiC(OOOl)), the interface is now well characterized. The 
first carbon layer, possibly with defects To*], strongly in- 
teracts with the substrate and forms covalent bonds. 
This layer yet acts as a buffer layer, so that the second 
graphitic plane exhibits graphene low-energy electronic 
properties 11 -17]. Thus, an overall linear dispersion is 
observed by angle resolved photoemission spectroscopy 
(ARPES), with a doping in electrons that shifts the Dirac 



point 0.4 eV below the Fermi energy [Tsi [l9|. 

Intentionally adsorbed species on the top of graphene 
can then be used to tailor its electronic properties. By de- 
positing molecules (20I [2T| . transition metals [22] or other 
elements[3 13: one can control the position of the 
Fermi level without altering the typical band structure of 
graphene. Adsorbates however are likely to act as scat- 
tering centers and increase the elastic scattering rate and 
then affect the transport properties [13, 113, [13 • For this 
reason, modifications of the interface seem more adapted 
since they occur further from the graphene layer. For 
example, gold can be introduced at the interface in or- 
der to dope the graphene layer with holes [1^. Recently, 
passivation of the substrate has been achieved using hy- 
drogen. In this case, the buffer layer is decoupled from 
the substrate and recovers a low-energy electronic struc- 
ture similar to pristine (neutral) graphene [26l [27| . 

On the C terminated face (SiC(OOOl)), the situa- 
tion is different. For ultra- high vacuum grown samples, 
coexisting native surface reconstructions of SiC (0001), 
namely the SiC(3 x 3)i2| and the SiC(2 x 2)c[2|, 
are found at the interface [s^, HJ. The graphene- 
substrate interaction is much weaker than on the Si 
face. Graphene is almost free-standing on the SiC (3 x 
3) reconstruction [s^] and weakly interacting with the 
SiC(2 X 2)c reconstruction [3^ 33]. Monolayer sam- 
ples show a band structure that resembles that of free- 
standing graphene, but with an intrinsic doping in 
electrons [34. The Dirac point is found 0.2 eV below the 
Fermi level^ [13] . For this system, effects of modifica- 
tions of the interface -either intentional or due to native 
defects- remain to be investigated. 

Here we focus on the SiC(2 x 2)c interface, which has 
a simple Si adatom and C restatom structure [2^. The 
electronic properties of graphene on SiC(2 x 2)c have al- 
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ready been investigated using ab initio calculations |33l|. 
For this system, without defects, we have shown that i) 
the total energy of the system is not sensitive to graphene 
versus SiC lattice translations, ii) the graphene layer lies 
far from the substrate (3.lA), iii)linear dispersion is pre- 
served in the vicinity of the Dirac point, iv) interaction 
with the adatom states occurs and results in a band an- 
ticrossing. It leads to band structure modifications for 
energies higher than 0.5 eV with respect to Ed and v) 
there is no charge transfer from the substrate. This lat- 
ter point is in apparent contradiction with experiments. 
Indeed, as mentioned above, transport measurements fl| 
and ARPES measurement s 1341 indicate a doping in elec- 
trons. In particular in ref. l34l . for a submonolayer sample 
with both SiC(3 x 3) and SiC(2 x 2)c types of interface, 
authors find the graphene Fermi level 0.2 eV above Ejj. 

In this paper, we show by means of ab initio calcula- 
tions that a class of interface defects- related to the spe- 
cific atomic structure of the interface- allows to change 
the electron doping level of graphene on SiC(2 x 2)c 
(G_2 X 2 in the following). These defects could be in- 
tentionally generated by chemical modifications but scan- 
ning tunneling microscopy (STM) data indicate that they 
are natively present on as-grown samples. Incidentally, 
defects observed by STM could explain the discrepancy 
between the calculated electronic structure for a perfect 
interface structure and the experimental results men- 
tioned above. 

Based on STM data, we propose in section II. C a 
structural model for the defects in which an adatom dan- 
gling bond is saturated by a H atom. We then address 
the effect of these defects on the graphene-substrate inter- 
action, and on the electronic properties of the graphene 
overlayer by means of ab initio calculations. Finally, we 
compare ab initio results to STM data in order to check 
the validity of the defect model. 



II. EXPERIMENT AND CALCULATION 
DETAILS 

A. Experimental aspects 

Substrates are n doped 6H-SiC(0001) (C face) polished 
by NovaSiC. Samples are prepared under ultrahigh vac- 
uum (UHV), following a procedure presented in ref. [13 
and [S^. Samples are first annealed at 850°C under Si 
flux in order to get a clean surface. Further annealing 
at 950 - 1000° C provides a weh ordered SiC (3 x 3). 
Finally, increasing the temperature to 1100° C leads to 
a graphene coverage of a fraction of a monolayer. The 
growth is controlled by low-energy electron diffraction 
(LEED) and Auger spectroscopy. For graphitized sam- 
ples, the LEED pattern shows two surface reconstruc- 
tions, the SiC (3x3) [Hi and the SiC(2 x 2)c[2i|. Already 
for a submonolayer coverage, the graphene signal is ring- 
shaped with modulated intensity, indicatiiig a str ong ro- 
tational disorder within the graphene film |30l432l |34|. 



Scanning tunneling microscopy (STM) measurements 
were performed in situ, under UHV at room tempera- 
ture. Tips are mechanically cut Ptir. STM gives ac- 
cess to the local density of states (LDOS) at the sur- 
face of the samples, with atomic resolution. In the case 
of graphene on SiC, it is well known that low-bias im- 
ages give access to the LDOS of the graphene layer while 
at high bias, graphene becomes transparent and images 
show the LDOS of the interface, i.e. under the graphene 
layer [TTi [35| . As presented in previous papers [lo, Is^ - 
the samples show bare SiC (3 x 3) domains and graphene 
monolayer islands on the SiC(3 x 3) and the SiC(2 x 2)c 
reconstructions. In agreement with LEED, there is a 
strong rotational disorder between the graphene islands 
which gives rise to various moire patterns (32|. 



B. calculational framework 

Calculations are carried out using the VASP codejs^, 
which is based on the density- functional theory (DFT). 
We use the generalized gradient approximation[37] to- 
gether with ultrasoft pseudopotentials (38j . The 4II-SiC 
substrate is modeled by a slab containing 4 SiC bilay- 
ers. The backside of the slab is passivated by hydro- 
gen (see Fig. [5] (b)). An empty space of 8 A separates 
the graphene layer from the next SiC slab. We use a 
plane wave basis cutoff of 211 eV. The ultrasoft pseu- 
dopotentials have been extensively tested[l^ [s^. Inte- 
gration over the Brillouin zone is carried out within the 
Monckhorst-Pack scheme, using a 6 x 6 x 5 grid. The K 
point at the corner of the Brillouin zone is included in the 
grid in order to get an accurate description of the band 
structure at the Dirac point. All the structures were fully 
converged, with residual forces smaller than 0.015 eV/A. 

STM images are simulated by cross sections of j^"!^ 
integrated over different ranges of energy, depending on 
the imaging bias. 



C. Modeling of the defects 

Fig. [1] (a) shows the ideal structure of the SiC(2 x 2)c 
reconstruction. A Si adatom sits in a three-fold coordi- 
nated hollow site (H3) and saturates 3 out of 4 C dangling 
bonds (DB) of the ideal surface[2'95. This leads to 2 DB 
per unit cell, one on the Si adatom and one on the fourth 
C atom (located at the corners of the diamond cell in 
Fig. 1 (a)). We call this fourth atom a restatom, by anal- 
ogy with the dimer adatom stacking fault (DAS) model 
that describes the Si(lll)(7 x 7)[39] surface reconstruc- 
tion, the SiC(2 x 2)c showing however a much simpler 
structure. From ab initio calculations 33j and STM mea- 
surements [13, [Si , it is established that a charge transfer 
occurs, resulting in an empty DB for the Si adatom and 
a filled DB for the C restatom. Such a charge transfer 
mechanism is also observed on the Si(lll)(7 x 7) and 
Ge(lll)c(2 X 8) surfaces, which involve restatoms and 
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adatoms[4ij. Fig. [U (b) and (c) are dual bias STM im- 
ages that illustrate this statement. In Fig. [T](b), a filled- 
state image, C restatoms appear at the corners of the 
dashed-line unit cell. In Fig. [T] (c) , an empty-state im- 
age, the Si adatom appears within the dashed-line unit 
cell. 

The real interface however exhibits a significant 
amount of defects. On Fig. [1] (d), a high-bias filled-state 
STM image of G_2 x 2, the SiC(2 x 2)c restatom lattice is 
interrupted by defects that look like "missing" restatoms 
surrounded by bright spots. These defects are still de- 
tected on low-bias images, together with the graphene 
lattice. In Fig. [T](e), the typical low-bias feature of a de- 
fect is seen on the left side of the image. Thus, the defect 
contributes to the DOS at the Fermi level {Ep), within 
the substrate surface bandgap, and could therefore be re- 
sponsible for the doping of the graphene layer. Note that 
the graphene island in Fig. [T] (e) presents the same ge- 
ometry as in the calculations (presented below), i.e. the 
graphene and the SiC surface lattices are aligned. Away 
from the defect, the graphene low-energy LDOS appears 
as a honeycomb lattice, modulated by a moire pattern 
which corresponds to the graphene 5x5 /SiC 4x4 com- 
mon cellfiol]. Its pseudo-cell is shown on the image. 

Although STM does not allow chemical characteriza- 
tion of the surface, it can provide useful hints on the 
structure of the defects in the real interface. Fig. [T] (b) 
and (c), show a typical defect, within the solid- line dia- 
mond cell. On filled-state STM images (Fig. [I](b)), the 
defects generally appear as a "missing" restatom, with 
one or more of its nearest adatom neighbors visible at 
both sample bias polarity (at the corners of the solid-line 
diamond cell). On empty-state STM images(see Fig. [T] 
(c)), we see the regular lattice of adatoms, with stronger 
signal on the 3 neighbors of the defect. These empty- 
state and filled-state features can arise from an impu- 
rity that saturates the restatom DB, and thus hinders 
the charge transfer from the neighboring adatoms. We 
propose a simple model to feature these defects that in- 
volves a H atom, which is also interesting since for the Si 
face of SiC, hydrogen decouples graphene from its sub- 
strate or saturates interface DB [42|. Note that such a 
cancellation of charge transfer has been observed on the 
Ge(lll)c(2 X 8) surface after hydrogen adsorption (is). 
In this case, the presence of the defect leads to a sim- 
ilar empty/filled states contrast on restatoms/adatoms, 
which further supports our model. 

In order to build the calculation cell, we start from the 
converged G_2 x 2 structure presented in ref [33] , namely 
a 5x5 graphene cell on top of a 4x4 SiC supercell, with- 
out rotation (see Fig. [2] (a)). The mismatch between the 
two supercells is smaller than 0.2 %. We then introduce 
some defects at the interface. One H atom is introduced 
in the G_2 x 2 structure, on top of the restatom, midway 
between the restatom and the graphene layer i.e. at 2 A 
from the restatom to be saturated. Since the supercell 
contains 4 SiC(2 x 2)c cells, there are 4 available loca- 
tions for the defect - labeled CI, C2, C3 and C4 in Fig. H 



(a). The 3 configurations CI, C3 and C4 appear similar 
with a defect located under a C-C bond. The defect in 
C2 is located under a graphene hexagon. All 4 configura- 
tions were studied and actually lead to similar electronic 
properties. We therefore present results only for the CI 
configuration, some results for the C2 configuration are 
given in the supplementary information . 




FIG. 1: (Color online) STM images of G_2 x 2 islands. Sample 
biases (Vs) are indicated on the images, (a) The SiC(2 x 2)c 
reconstruction. The unit cell is shown in dashed lines. The 
size of the atoms refers to their height in the cell. Only the 
topmost SiC bilayer and the Si adatoms are represented, (b) 
and (c) 4x4 nm^ dual high-bias STM images. The SiC(2x2)c 
appears. The undisturbed unit cell is given in dashed lines. A 
typical interface defect is shown at the center, highlighted by 
the full line cell, (d) 20 x 20 nm^ high-bias STM image. Point 
defects are seen at the interface, (e) 5 x 5 nm^ low-bias STM 
image of a G_2 x 2 island with no rotation with respect to the 
substrate. A defect is seen on the left side of the image. Away 
from the defect, at the atomic scale, the graphene honeycomb 
lattice is seen. At a larger scale, a periodic pattern (5x5 with 
respect to the graphene lattice) is visible. Its pseudo-cell is 
shown in full line. 



III. 



RESULTS 



A. Atomic structure 

After relaxation, the H atom has indeed formed a bond 
with the restatom and is now located only 1.10 A above 
it. Some other modifications of the atomic positions are 
also noticeable, but they do not significantly depend on 
the defect location. The mean graphene layer-adatom 
distance (initially of 3.10 A) dropped by 0.15 A. On the 
lower half of the 5x5 graphene supercell, atoms of the 
graphene layer have gone down toward the substrate, 
slightly increasing the already existing layer corrugation, 
from 0.15 A peak to peak for the defect-free structure to 
0.19 A. Atomic displacements are thus small and range 
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between 0.05 A and 0.13 A. Maximum displacements 
are observed for graphene atoms near the adatoms first 
neighbors of the defect. Conversely, adatoms of the re- 
construction have rose toward graphene, with a maxi- 
mum displacement of 0.11 A for the 3 adatoms neigh- 
boring the defect (for the CI configuration: the adatoms 
labeled 1, 3 and 4 on Fig. [H (a)). The 4th adatom ex- 
hibits a displacement of 0.06 to 0.08 A, depending on the 
location of the defect. 

Regarding the stability of the structures, the location 
of the defect does not seem to be a decisive factor. The 
total energy of CI, C3 and C4 differ only by maximum 
2 meV and C2 differs from the 3 other configurations by 
maximum 9 meV per unit cell. The less stable configu- 
ration is C2. 

Thus, the presence of defects implies some structural 
modifications. The graphene-substrate distance tends to 
decrease, which could come with modifications of the 
electronic properties of the graphene layer and especially 
a strengthening of the graphene-substrate interaction, 
contrarily to what is observed on the Si face samples 
when H is introduced at the interface. We will address 
this question in the following. 



B. Electronic structure 

In the previous study [111, for the bare SiC(2 x 2)c 
surface in the 4x4 geometry, the integrated DOS showed 
two narrow sets of bands corresponding to the adatom 
and restatom states, with the Fermi level at the top of 
the restatom band. Adding a graphene layer on top of 
the reconstruction did not qualitatively change the DOS 
of the interface (see the DOS of the defect free G_2 x 2 
structure in Ref. [H, Fig. SI (b)). When a defect is 
added to the structure - regardless of its location -, the 
total DOS shows a marked shift of the bands toward filled 
states compared to the DOS of the bare SiC(2 x 2)c or 
the G_2 X 2(see Fig. [2] (c) for the CI configuration and 
Ref. 13, Fig. SI (a) for the other configurations). 

A detailed observation of the band structure provides 
more information on this aspect. The band structure of 
the CI system is represented in Fig. [2] (d). First of all, 
the TT bands of graphene are still present. The Dirac point 
is located approximately 0.48 eV below the Fermi level, 
indicating a significant electron doping. Bands are linear 
over +/—100 meV with respect to the Dirac energy (Ed). 
For higher energies with respect to Ed, a band anticross- 
ing effect appears, involving the tt* band and an adatom 
band, which distorts the graphene band around Ep. 
Thus, a hybridization of graphene states with adatom 
DBs occurs. This feature was already present in the de- 
fect free system but here, the energy splitting at anti- 
crossing increases from AE « 0.35 eV (defect-free struc- 
ture) to AE « 0.45 eV (see Ref. H Fig. S2(a)) confirm- 
ing an increased interaction with graphene. Regarding 
the TT band, interaction with states from the C restatoms 
remains very weak, the restatoms still being located far 



from the graphene layer. Finally, very small or vanishing 
gaps are observed at the Dirac point. CI shows a gap 
of « 15 meV (see Fig. 2 (d)) and the maximum gap of 
« 30 meV is obtained for C2 (see Ref. Fig. S2(b)). 
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FIG. 2: (Color online) A 5 x 5 graphene cell on top of a 4 x 4 
S'iC(OOOl) substrate in the (2 x 2)c reconstruction with one 
defect (H atom) per cell, (a) Top view of the supercell. Only 
the last SiC bilayer, the Si adatoms, the H atom and the 
graphene atoms are represented. Within the supercell (black 
diamond cell), all atoms are plotted once and their size refers 
to their height. The Si adatoms are numbered from 1 to 4 and 
the CI to C4 labels refer to the 4 possible locations of the H 
atom, (b) Side view of the supercell after relaxation for the CI 
configuration, (c) density of states for the CI configuration, 
(d) Band structure for the CI configuration in dotted line. 
The dispersion of an isolated graphene layer, shifted to align 
Eds, is given in solid line. The inset shows a zoom of the CI 
configuration band structure in the vicinity of the K point. 



C. Partial charge distribution versus STM images 

In order to check if our defect model accurately de- 
scribes the real system, it is interesting to compare par- 
tial charge density maps to STM data. We however stress 
that such a comparison is limited to qualitative consider- 
ations because of the following reasons. The real doping 
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is certainly smaller than the one we get from ab initio 
calculations, due to a smaller density of defects in the 
real system. Thus, the relative positions of Ejj and Ep 
are not quantitatively reproduced. Moreover, it is known 
that the DFT gap is systematically underestimated and 
the position of the dangling bond states bands is therefore 
only qualitative. Furthermore, the extension in vacuum 
of the numerical wavefunction is significantly underesti- 
mated. Thus, the low-bias STM images are compared to 
cross sections taken just above the graphene atoms. For 
the same reason, high-bias STM images - which reveal 
the SiC surface states - are compared to cross sections 
taken just above the adatoms, i.e. below the graphene 
layer. Note finally that, due to the limited size of the 
supercell considered here, the calculation essentially pro- 
vides information on the vicinity of the defect. Away 
from the defect, the STM image is reproduced by the 
defect-free interface modelfsS*]. 

Fig. 3 (a) is a recall of the calculated structure (trans- 
lated with respect to Fig. 2 (a) to get the defect in 
the middle of the cell). In the first instance, we concen- 
trate on the states belonging to the substrate, which are 
probed on high-bias STM images. Thus, partial charge 
density maps are taken just above the adatoms (see the 
right part of Fig. 3 (a)). Fig. 3 (b) is a cross section 
of l^'p above the adatoms, integrated from —1.0 eV to 
0.0 eV. It has to be compared to high- bias fiUed-state 
STM images like in Fig. [T] (b). As expected for the 
SiC (2 X 2)c reconstruction, states are present on the re- 
statoms. However, at variance with the ideal reconstruc- 
tion, a restatom is "missing" (indicated by the arrow) . It 
corresponds to the one bonded to the H atom, whose 
DB has been saturated by the H atom. Moreover, states 
are also detected on the adatoms, with higher intensity 
-even brighter than on the restatoms- on the 3 neighbors 
of the defect, labeled 1, 3 and 4 on Fig. 3. This higher 
intensity is in good agreement with the STM data in the 
vicinity of the defect site, except that calculations show 
similar charge density on the 3 neighbor adatoms while 
STM images usually do not. We will comment on this 
later. 

Fig. 3 (c) also shows a cross section above the adatoms, 
but for l^fp integrated from 0.0 cV to +0.5 eV, which 
has to be compared to high-bias empty-state images like 
in Fig. [1] (c). In this case, states are localized on the 
adatoms, with equivalent intensity. In the empty-state 
STM image of Fig. [1] (c), we also have signal only on 
adatoms, but with increased intensity on the neighbors 
of the defect. We however observed on positive sample 
bias STM images that the contrast near defects depends 
on the probed energy window. The larger the energy 
window, the less noticeable is the perturbation from the 
defect. 

We now focus on the states belonging to the graphene 
layer. Thus, as illustrated in Fig. 4 (a), we consider par- 
tial charge density maps taken on a plane located just 
above the graphene atoms (Fig. 4 (b) and (c)), and also 
above the adatoms when the graphene-substrate interac- 
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FIG. 3: (Color online) Partial charge density maps for the 
CI configuration displaying the substrate states (Contrast is 
reversed with respect to Fig. 1: higher density is darker), 
(a) (left) Calculated CI structure. The cell is translated with 
respect to Fig. 2 (a) so that the H atom, highlighted by an 
arrow, is in the middle of the image. The size of the atoms 
refers to their height in the supercell. The Si adatoms are 
numbered from 1 to 4, adatoms 1, 3 and 4 being the nearest 
neighbors of the defect. (Right) Side-view of the calculated 
structure, the dotted line indicates the section plane of the 
partial charge density maps, (b) and (c) Maps of l^*]^ inte- 
grated from -1.0 eV to 0.0 eV and from 0.0 eV to -1-0.5 eV 
resp. 



tion is considered (Fig 4 (d)). In the following, we will see 
how features characteristic of the defect free G_2 x 2 (see 
Fig. 1 (e), away from the defect) are reinforced due to an 
increase of the interaction, that is: supercell dark/bright 
halves. 3- fold symmetric line pattern and switched-off 
atoms j33|. 

Fig. 4 (b) shows a cross section of l^'P, just above 
the graphene layer, integrated from —0.55 eV to —0.2 
eV. Notice that, at variance with Fig. 3, four supercells 
are represented. A schematic view of the data is given 
on the right part of the image, in which we also added 
the Si adatoms and the H defects as guides to the eye. 
As expected for the considered energy range, i.e. where 
the graphene/substrate interaction remains weak (energy 
window in between the SiC (2 x 2) dangling bonds states. 
Fig. 2 (d)), we observe essentially a honeycomb con- 
trast. Within the supercell, which is indicated by a black 
diamond cell on Fig. 4 (b), the upper left half of the su- 
percell appears higher (darker with this color scale) than 
the lower right one. This feature arises from the topo- 
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graphic modulation described earlier in this paper. If we 
then consider a larger area than the supercell, we notice 
that the higher graphene atoms form a 3-fold symmet- 
ric pattern made of lines of hexagons (see the schematic 
view on the right part of Fig. 4 (b)). 

In Fig. 4 (c) and (d), cross sections of integrated 
from —0.1 eV to -1-0.1 eV are represented. In this en- 
ergy range, the coupling between the electronic states of 
the graphene layer and the substrate states, namely the 
adatoms states, is strong (discussed in section III. B, Fig. 
[2] (d)). The left part of Fig. 4 (c) shows data taken just 
above the graphene layer, illustrated by a schematic view 
on the right part. The supercell is again indicated by a 
black diamond cell. In the schematic view, Si adatoms 
and H defects are added to the data. The graphene atoms 
within the supercell are represented by small black balls 
when they display significant charge density and by dot- 
ted circles when the charge density is much smaller. In 
the lower part of the 5x5 graphene supercell, we clearly 
see a contrast asymmetry between the graphene A and B 
sublattices. It seems to arise from the symmetry of the 
local stacking. Indeed, in this region, 3 graphene atoms 
belonging to the same graphene sublattice are located 
near adatoms and are switched off. In the upper part 
of the cell, one graphene atom, belonging to the other 
graphene sublattice and located on top of an adatom is 
also switched off. 

If we then move down to the top of the adatom with 
Fig. 4 (d), charge density appears on the adatoms, 
mainly on the 3 neighbors of the defect. Actually, Fig. 4 

(c) and (d) give a real space image of the interaction be- 
tween graphene and the reconstruction adatoms. Fig. 4 

(d) shows that hybridization occurs mainly between the 
TT* band of graphene and the DBs of the 3 adatoms neigh- 
boring the defect. Finally, notice that high-bias empty- 
state images like in Fig. [1] (c) in fact correspond to an 
average of data in Fig. 3 (c) and 4 (d), with a DOS signal 
on the Si adatoms, which is increased on the 3 neighbors 
of the defect. This illustrates the fact, evoked in the dis- 
cussion in connection with Fig. 3 (c) that the aspect of 
the adatoms depends on the integration window. 



IV. DISCUSSION AND CONCLUSION 

To sum up, for the kind of defect we study, a restatom 
DB is saturated by an H atom. As a result, the charge 
transfer from the 3 first neighbors adatoms toward the 
restatom is impossible, which leaves a delocalized elec- 
tron on these 3 adatoms. This local perturbation leads 
to a doping in electrons of the graphene layer. The Dirac 
point is namely shifted by « 0.5 eV under the Fermi 
level, which would correspond to n ~ 2 .lO^'^cm"^ for 
ideal graphene. Such a doping is of the order of the elec- 
tron density created by the defect(rt « 7.6 .IO^'^^cto"^ i.e. 
one electron per defect). Modifications of the electronic 
structure come with a maximum reduction by 0.2 A of the 
graphene- adatom distance. The graphene-substrate in- 



b and c 
-a- --«-d 




-0.1/+0.1 eV 



FIG. 4: (Color online) Partial charge density maps for the 
CI configuration showing the substrate states (Contrast is 
reversed with respect to Fig. 1: higher density is darker). 
States from the graphene layer are considered here, (a) Side- 
view of the calculated structure, the dotted-lines indicate the 
height of the resp. partial charge density maps showed in 
the following, (b) (left) Map of I'I'p integrated from —0.55 
eV to —0.2 eV (low-interaction region), taken just above the 
graphene atoms. The black diamond cell indicates the 5x5 
graphene supercell, thus, 4 of them are represented on the 
image, (right) Corresponding schematic view. The honey- 
comb lattice represents the graphene lattice, higher density 
is darker. The Si adatoms and the H defects are also repre- 
sented, (c) and (d) Maps of j^tp integrated from —0.1 eV to 
-1-0.1 eV (interacting region) taken just above the graphene 
atoms and just above the adatoms resp. In (c) the data (left) 
is also completed by a schematic view (right). Within the 
supercell, the C atoms belonging to the graphene plane are 
depicted as small black balls and as dotted circles (see text). 
Again, Si adatoms and H defects are represented. In (d), the 
position of the H defect is highlighted by an arrow. 



teraction is slightly enhanced, and involves the 3 adatoms 
neighboring the defect. However, the linear dispersion is 
preserved for energies within +/- lOOmeV with respect 
to the Dirac point. 

Based on STM results, our model seems to repro- 
duce well the perturbations induced by the defect to the 
SiC(2 X 2) . Calculations however show defects with higher 
symmetry than in the real system since the 3 adatoms 
neighboring the defect show similar charge density(Fig.3 
(b)). Another calculation with a shift in the relative po- 
sition between graphene and SiC lattices, while leading 
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to a qualitatively similar electronic structure, results in 
a less symmetric defect in the direct space (see Ref. |3, 
Fig. S3). Moreover, we actually considered the most sim- 
ple defect configuration involving an H adsorbate while 
in the real system other species might come into play. 
Additionally, the system studied exhibits a high density 
of defects (7.6 lO-'^^cm"^) due to the technical constraint 
that limits the size of the supercell. Such a high defect 
density might impose symmetries in the system. 

Considering the graphene-substrate interaction, be- 
cause of the high density of defects in the calculated 
structure, the doping and the enhancement of the in- 
teraction with the substrate are most probably overes- 
timated with respect to the real system. Other groups 
have measured an intrinsic doping in electrons of n « 
10^^cm~^[9, 34]. From the typical STM image in Fig. 
[T] (a), we find a defect density of n w lO^^cm^^. The 
experimental defect density is high enough to generate 
the experimentally measured electron concentration and 
defects are thus a plausible cause for the doping of the 
graphene layer. Note that on STM low-bias images, no 
(V3 X V3)i?30° standing wave pattern are observed in 
the vicinity of defects. Thus, the defects -located at 3A 
under the graphene layer- do not lead to intervalley scat- 
tering. 

Finally, in our system, modification of the interface 
by introducing H atoms {10^*cm~^) leads to a doping 
in electrons of the graphene layer and also slightly en- 
hances the graphene-substrate interaction. Riedl et al. 
in Ref. [1^ have shown that, in the case of graphene 



grown on the Si face of SiC, the graphene-substrate in- 
terface can be H passivated and the graphene-substrate 
interaction vanishes. Starting from the strongly cou- 
pled buffer layer, after H treatment, they obtain a quasi- 
free-standing graphene layer, neutral, with linear disper- 
sion. Thus, the molecular H migrates under the buffer 
layer and bonds to Si atoms of the substrate, break- 
ing the covalent bonds between the buffer layer and the 
substrate [27|. As a result, graphene characteristics are 
restored in the first graphitic layer. In the case of the 
C face, we conclude from our results that in order to 
decrease the graphene-substrate interaction, the density 
of H atoms introduced has to be sufficient to destroy 
the SiC(2 x 2)c reconstruction and saturate the ideal 
SiC surface. Otherwise the presence of H atoms would 
rather lead to a doped layer with an enhanced graphene- 
substrate interaction. 
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